Heating the O(N) nonlinear sigma model* 



Harmen J. Warringa 

Department of Physics and Astronomy, Vrije Universiteit, 
De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 

February 1, 2008 



Abstract 

The thermodynamics of the 0{N) nonlinear sigma model in 1 + 1 dimensions is studied. We 
calculate the finite temperature effective potential in leading order in the l/N expansion and 
show that at this order the effective potential can be made finite by temperature independent 
renormalization. We will show that this is not longer possible at next-to-leading order in l/N. 
In that case one can only renormalize the minimum of the effective potential in a temperature 
independent way, which gives us finite physical quantities like the pressure. 
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1 Introduction 



The nonlinear sigma model is a scalar field theory with an O(N) symmetry. It is described by a 
Lagrangian density which only consists of a kinetic term, 

£ = \d^df<k , (i) 

and a constraint which enforces all the <j) fields to lie on a N — 1 sphere: 

Ux)(t>i{x)=N/g 2 i = l...N. (2) 

This model has some nice features in 1 + 1 dimensions, which makes it interesting to study as a 
toy model for QCD. First it is renormalizable. Furthermore it is asymptotically free, such that 
at very high temperatures it approaches a free field theory. The model also has a dynamically 
generated mass for the (ft fields. If N = 3 the model has instanton solutions. Finally, for N = 2 
we recover a free field theory, which can be used as a check of the calculations. 

In this article we will study the thermodynamical properties of the nonlinear sigma model. In 
particular we will calculate the pressure. In Sec.|2]we briefly discuss some aspects of thermal field 
theory. In Sec. 01 we calculate the pressure in the weak-coupling expansion. In Sec.@J we calculate 
the effective potential and pressure to leading order in the 1/N expansion. The next-to-leading 
order (NLO) correction is discussed in Sec. |5] 



2 The pressure in a field theory 



In this section we briefly review how one calculates the pressure in a thermal field theory. For a 
more complete introduction see Refs. pJEj. 

In classical statistical mechanics one can derive all thermodynamic quantities from the parti- 
tion function. The partition function Z is given by 



(3) 



where the sum is over all eigenstates of the Hamiltonian H and (3 = 1/T, the inverse temperature. 
For example the pressure V is given by 



= 1 d log Z 

(3 dV 



(4) 



We next express the partition function in terms of fields. The easiest way to do this is to consider 
a transition matrix element in ordinary field theory. One can write such a transition element in 
terms of a path integral in the following way 



exp[-i(tf - ti)H] 



T>4> exp 



dt / d d xC( 



(5) 



where £ is a Lagrangian density which has a Minkowskian metric and does not have derivative 
interactions. Now if one makes the identification t = —it one finds 



0/ ex V [-PH] 



Vcj) exp 



dr / d a x £((/)) 



(6) 
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where we from now on denote the zero component of a (d + l)-vector by r and hence use a 
Euclidean metric. The last equation enables us to write the partition function in terms of a path 
integral, 



= j 'vcj) exp - ^ 'dr J d d x £(0) 



J ^( t =0)=</»(t=/3) 



(7) 



where one implicitly integrates over all states which obey the periodicity condition 0(r = 0, x) = 
4>(t = f3,x). So equilibrium thermal field theory is in essence a Euclidean field theory where 
one dimension (r) is compactified to a circle. As a consequence, the Fourier transform of a field 
becomes a sum over modes, 



1 



d d k 
(2vr) d 



iuj n T+ik-x 



T+ik-x 



(8) 



where uj n = 2imT. This implies that in a loop diagram one should not take the integral over 
internal momentum but rather the sum-integral ^. 

Now for example the partition function of the nonlinear sigma model is given by 



r N r r 13 r 

Z = [Jp^^^W-JV/^exp - dr dxC 

i=l x I JO J 



(9) 



where from now on we work in one spatial dimension. To obtain the pressure we have to calculate 
Z. We will follow two paths. The first one is making an expansion around g 2 = 0. This will 
only give us the leading term of the pressure. The second way is an expansion in 1/N which will 
generates additional contributions which are non-analytical in g 2 . 



3 The pressure in the weak-coupling expansion 

One can get rid of the constraint by integrating out one of the (j) fields, which results in 

.N-l p .p . 

Z= H Vm fj 6 (N/g 2 - TTiiTi) exp - dr dx C cS (tt) 

J i=l x I JO J 



(10) 



where 9(x) is the step function and the effective Lagrangian density C e g is given by 

1. ... 9 2 (ir^TTj) 2 1 



£eff(7r) = + y N-l 2 ^ ~ 2 PVlog ( N/g2 ~ ^ 



(11) 



For small values of g 2 the 9(x) function is only vanishing when n(x) is large. Since large values of 
7r give a small contribution to the path integral we approximate 9(N/g 2 — TTiTTi) « 1 which gives 



Z 



.N-l r .0 . 

J IJ Vmexp^- dr J d d x£ cS (TT) 



(12) 



We will not calculate Z but rather -^y log Z, where V is the volume of our 1 dimensional space. 
Because logZ is an extensive quantity, i.e. it is linear in V, the pressure is equal to -^y\o<gZ. 
Since in general -gy log Z does not vanish at zero temperature, we subtract the zero temperature 
contribution to normalize the pressure to zero at zero temperature. 
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If g 2 = it can be seen from C e $ that one has N — 1 noninteracting tt fields. Hence it is easy 
to show that leading term is equal to the pressure of N — 1 free fields 



V 



N — 1 



K 



\og{K 2 



K 



\og(K 2 



(N 



1)*-T 2 



where K 



(io n , k ) is a Euclidean two vector and we defined 

d 2 k 
(2vr) 2 • 



A' 



(13) 



(14) 



By calculating the loop diagrams, one can show that up to and including order g 4 one only 
finds the pressure of a free gas in <i = 1 + 1 [HIE]- However one finds corrections to the free 
pressure in a 1/N expansion. This may indicate that the pressure is completely non-analytical 



m g 2 . 



4 The effective potential in leading order in 1/N 

Another way to implement the constraint on the (ft fields is by using a Lagrange multiplier field 
which we will denote by a. This gives the following expression for the partition function, 



Z = j JJPtfoPaexp j-i j dr J dxd^ 



P 

dr I dxa(x)[Mx)Mx)- N/g 2 }} . (15) 



o 



In this way the action still is quadratic in the (ft fields, so one can easily integrate them out. This 
gives 



Z = J Vaexp / — S[a] + ^-^ J dr J dxa(x)\ 



(16) 



where 

S[a] = y Tr log[-<9 2 + a(x)\ . (17) 

The pressure is equal to the minimum of the effective potential, which one can calculate by 
expanding the a field around its vacuum expectation value m 2 . By considering the propagator 
of the (ft fields, one can show that to leading order in 1/N, m is equal to the physical mass of 
the (ft fields. This is, however, not longer the case at NLO, [51 . The effective potential can be 
obtained from the effective action by division by (3V . To calculate the effective potential we write 
a = m 2 + aj and expand the action around m 2 [Zj , 

S[a] = ^Tr log[-a 2 + m 2 ] + ^Tr (3^2 *) + \^ «) ' + 0(1/ V#) ■ (18) 

From this equation it can easily be seen that the effective potential can be calculated in a 1/N 
expansion. The leading order effective potential is given by the classical action. The corrections 
are obtained by integrating over the a field. 

To calculate the leading order effective potential we introduce a momentum cutoff A and 
subtract m and T-independent constants from the effective potential. This subtraction will not 
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change the physics, since it only shifts the whole effective potential by a constant. One finds for 
the effective potential at leading order in 1/N 

N [jf p log(P 2 + m")-J p lo g (P 2 

,2 / \ rp2 



V(m 2 



N 
~2 



N 
~2 



m 

n 



rn 

47T 



1 + log • 



m z 



+ — J ((3m) 



(19) 
(20) 



where g\, is the bare coupling constant. Jo (/3m) is given by 

p 2 n(u p ) 



J (/3m) 



8 



dp '■ 



(21) 



where n(x) = l/(e x — 1) and to 2 = p 2 + m 2 . One is able to renormalize the leading order effective 
potential in a temperature independent way by replacing g 2 — > Z g 2g 2 (fi) where 



Z g 2 



1 + f- log - 



and g 2 = g 2 {p)- From this equation it follows that the /3-function of g 2 is given by 



(3(g 2 



,dg 2 M_ g 4 (/i) 
d/i 2-7T 



(22) 



(23) 



The leading order /3-function is exact in g 2 . Since the /?- function is negative, g 2 approaches zero 
for large values of /i. This shows that the theory is asymptotically free. 

With use of the renormalization of the coupling constant one finds the following finite expres- 
sion for the effective potential 



V(m 2 ) = - 



1 



2 2/ 2 

m m / , a . ...... 



(24) 



One can easily show that the effective potential is independent of the renormalization scale \i. 
This is expected since the choice of \i is completely arbitrary. 

To obtain the pressure, one has to minimize the effective potential with respect to m 2 . Mini- 
mization gives the so-called gap equation 

,2 \ i i / ,,2 

^2 



1 



r 



i 



- „ > , 2 = J" lQ g ( ^2 J + J" J l (^ m ) = J" l0 § ( ^? 

P ^ + m z 47r \m z J 47r 47r \m 



(25) 



where Ji(/3m) is defined by 



Ji(/3m) 



dp 



n(u} p ) 



(26) 



The solution of the gap equation determines the leading order physical mass of the 4> fields as a 
function of temperature. At T = one can solve this equation to show that the mass is completely 
non-analytical in g 2 , 

m T=0 = nexp (^~~T^j • ( 27 ) 

We can use Eq. (|27jl to normalize the minimum of the effective potential at T = to be zero 
which gives 



V(m 2 



N 
~2 



m 



m 
Air 



1 + log 



(28) 
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The effective potential as a function of m for different temperatures is shown in Fig. (|T|). for the 
arbitrary choice g 2 (n = 500) = 10. The quantities T, in, fi, V/T and V/T are all in the same 
arbitrary units. The solid curve which is the minimum of the effective potential is equal to the 
pressure. 
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Figure 1: The leading order effective potential as function of m for different temperatures with 
g 2 (u = 500) = 10. 



5 Next-to-leading order correction in 1/N 

The term linear in a in Eq. ()18j) gives no contribution to the effective potential since it gives rise 
to a tadpole p]. The first 1/N correction to the effective potential stems from the last term of 
Eq. (j!8|) . By going to momentum space one can show that the correction is given by 



Vi(m 2 ) 



WMt 



i 



i 



2if P "° [% Q 2 + m 2 {P + Q) 2 + m 2 
We calculated this correction in Ref. In the limit A — > oo, one obtains 



(29) 



Vi(m 2 



-— (A 2 In In ^ 

8ir V rn z 



A 2 

m li —z 



m" /, . A 2 
In In — 



-in 



In 



A 2 \ 



4m 2 



) + F(m,T) , (30) 



where fh is defined in Eq. ()25|) . In (|30() . we have subtracted m and T-independent constants and 
dropped terms that vanish in the limit A — > oo. F(m, T) is a finite term and the logarithmic 
integral li is defined by 



P I dt-^— 
lo log* 



(31) 



where V stands for principal value. The first two terms of Eq. (|30|) are problematic. It is 
impossible to remove these divergences by renormalizing g 2 in a temperature independent way 
or by subtracting m and T-independent constants. However this is possible at the minimum of 
the effective potential. At the minimum, one can use the leading order gap equation, Eq. (|25|) . 
to show that fh is independent of T. So one could add 

A 2 f 4vr / 4tt\ , [ /4tt\1 1 , . 

In — o — exp T h exp > (32) 

9b \ 9 2 J L V5 2 yJJ 



8^ 
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to the effective potential which yields an effective potential that can be renormalized at the 
minimum. Using this renormalization at the minimum we have calculated the pressure V as a 
function of N. The result is depicted in Fig. ([%]). One clearly sees a crossover which is not a phase 
transition. This is in accordance with the Mermin-Wagner-Coleman theorem [SJ^j which forbids 
spontaneous breakdown of a continuous symmetry in 1 + 1 dimensions. The figure furthermore 
shows that the 1/N expansion is relatively good since the corrections are really of order 1/N. 
Finally it can be seen from the figure that the theory is asymptotically free, because in the limit 
T — > oo the pressure approaches the pressure of a free gas, Eq. I|13j). 
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Figure 2: Pressure V normalized to NT 2 as a function of temperature for different values of N 
with g 2 (n = 500) = 10 0. 



6 Summary and Conclusions 

We find that the pressure of the nonlinear sigma model in the weak-coupling expansion through 
order g 4 only consist of the free term. Furthermore, we showed that in a 1/N expansion we 
can renormalize the leading order effective potential in a temperature independent way. This is, 
however, impossible for the effective potential at next-to-leading order in 1/N. In that case one 
can only renormalize in a temperature-independent way a physical quantity, like the pressure. 
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